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Abstract 

The magnetic monopole condensate is calculated in the dual Monopole Nambu- 
Jona-Lasinio model with dual Dirac strings suggested in Refs.[l,2] as a functional of 
the dual Dirac string shape. The calculation is carried out in the tree approximation 
in the scalar monopole-antimonopole collective excitation field. The integration over 
quantum fluctuations of the dual-vector monopole-antimonopole collective excita- 
tion field around the Abrikosov flux line and string shape fluctuations are performed 
explicitly. We claim that there are important contributions of quantum and string 
shape fluctuations to the magnetic monopole condensate. 
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1 Introduction 



In Refs.[l,2] there has been suggested the dual Monopole Nambu-Jona-Lasinio (MNJL) 
model with dual Dirac strings as a continuum analogy of Compact Quantum Electrody- 
namics (CQED) which is defined for lattices as nonlinear U (1) gauge theory. It has a 
confining phase like QCD [3] and realizes confinement of "color" electric charges. Thereby, 
the investigation of CQED should help us to understand quark confinement. As has been 
shown in Ref.[4] the nonperturbative vacuum of CQED behaves like an effective dual 
superconductor with magnetic monopoles. Due to magnetic monopoles the electric flux 
between quarks rearranges and looks like the field produced by a dual Dirac string. As a 
result quarks interact via a linearly rising potential [5,6] that realizes quark confinement 
[7,8] and spontaneous breaking of chiral symmetry (SB%S) [8]. 

The NJL model [9] can be regarded as some kind of relativistic extension of the 
BCS (Bardeen-Cooper-Schrieffer) theory of superconductivity [10]. It also possesses a 
nonperturbative vacuum with a ground state of the same kind as in a superconductor 
in the superconducting phase. The latter has been the promoting idea of Refs.[l,2] to 
put the NJL-model to the foundation of a continuum space-time model realizing non- 
perturbative phenomena of CQED. 

The MNJL-model is based on the Lagrangian, invariant under "color" magnetic U(l) 
symmetry, of massless magnetic monopoles, self-coupled through local four-monopole 
interaction [1,2]: 

C(x) = x{x)il»d, X {x) + G[ X (x)x(x)} 2 - G 1 [x(x) 1 ,x( x )}\x( x h fl x( x )}, (1.1) 

where x( x ) i s a massless magnetic "color" monopole field, G and G\ are positive phe- 
nomenological constants responsible for the magnetic monopole condensation and the 
dual- "color" vector field mass, respectively. 

The magnetic monopole condensation accompanies the creation of massive magnetic 
monopoles Xm(x) with mass M, xX^ co ^ ec ^ ve excitations with quantum numbers of scalar 
Higgs meson field a with mass M a = 2 M and a massive dual-vector field C M with mass 
M c defined as [1,2]: 

M c = £;- £*IUM) + M2 HM% (1.2) 

where J\{M) and Ji{M) are quadratically and logarithmically divergent integrals [1,2] 
J '< M > = /^M^F= A2 - M2fa ( 1 + ^)' 

hiM) = / ^ijM^Wf = 4 1 + w) ~ (L3) 

Here A is the ultra-violet cut-off. The mass of the massive magnetic monopole field 
Xm(x) obeys the gap-equation [1,2]: 

M = -2G < x(0) X (0) >= f^-/i(M). (1.4) 

After the integration over magnetic monopole degrees of freedom the effective Lagrangian 
containing quarks, antiquarks and the fields of scalar a and dual-vector C M collective 
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excitations reads 



C cS (x) = ±F llv (x)F> u '(x) + ~M*C ll {x)C> t {x) 



+ l -d,o{x)d»a{x) - \mIo\x) 



a(x) 1 " 

K- 



M„ 



~\~ £free quark(^)) (l - ^) 

The coupling constants g and k are related by the constraint 

2 2 

' ;J 2 {M) = J 2 (M) = 1 (1.6) 



12vr 2 v ' 8tt 2 



or = 2^ 2 /3 [1,2]. 

£free quark (x) is the kinetic term for the quark and antiquark, 

Aree quark(x) = - £ TXl, J dr ^J^ ^^ 9^ ' *«(x - X 4 (t)). (1.7) 

We consider quark and antiquark as classical point-like particles with masses m q = 
rriq = m, electric charges Q q = —Qq = Q, and trajectories X"{r) and X?(t), re- 
spectively. The field strength F^(x) is defined [1,2] as F^(x) = £^(x) - *dC^(x), 
where dC^(x) = d»C u {x) - d u C^(x), and *dC v {x) is the dual version, i.e., *dC^(x) = 
\e^dC aP {x) (e 0123 = 1). 

The "color" electric field strength £^ v (x), induced by a dual Dirac string, is defined 
following [1,2] as 

where X 11 = X M (r, er) represents the position of a point on the world sheet swept by 
the string. The sheet is parametrized by the internal coordinates — oo < r < oo and 
< cr < 7r, so that X^(t, 0) = X^q(t) and X M (r, 7r) = Xq(t) represent the world lines of 



an antiquark and a quark [1,2,5]. Within the definition Eq. fll.8Q the tensor field 8 tll '{x) 
satisfies identically the equation of motion, d^F^ix) = J v {x). The electric quark current 
J u (x) is defined as 

J v {x) =Y.Qij dT^^-S^ix-Xiir)). (1.9) 

i=q,q T 

Hence, the inclusion of a dual Dirac string in terms of E^ v [x) defined by Eq. fll.8j ) satisfies 
completely the electric Gauss law of Dirac's extension of Maxwell ' s electrodynamics. 

As has been shown in Refs.[l,2] the vacuum expectation values of time-ordered prod- 
ucts of densities expressed in terms of the massless-monopole field, i.e., the magnetic 
monopole Green function 

G (Xt, ...,X n )=< 0|T(x(xi)rix(xi) . . . X( x n) T„X (^n))|0 >conn. , (1-10) 
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where r,(z = 1, . . . , n) are the Dirac matrices, are given by [1,2] 



G{xi, ...,x n )=< 0|T(x(a;i)rix(a;i) . . . x(x n )T n x(x n ))\0 > c 
= (M) < Q\t(xm(xi)TiXm(30i) . . . XM{x n )T n XM(x n ) 



i / d A x{-gxM{x)l v XM{x)C v {x) - kxm{x)xm{x)u{x) + C int [a(x)]}, (M) 
x e J |0 >\ZL ■ (1.11) 



Here |0 >( M ) is the wave-function of the nonperturbative vacuum of the MNJL-model 
in the condensed phase and |0 > the wave-function of the noncondensed perturbative 
vacuum. £ int [cr(a;)] describes self-interactions of the cr-field: 

C- mt [a(x)\ = -KM a a\x) - ^cr 4 (x). (1.12) 

The self-interactions £ in t[cr(x)] provide cr-field loop contributions and can be dropped out 
in the tree cr-field approximation accepted in Refs. [1,2]. In the tree cr-field approximation 
the r.h.s. of Eq. ( 1.1 If ) acquires the form 

G(x h . . . , x n ) =< 0|T(x(xi)rix(xi) . . . X(x n )T n x(x n ))\0 > coim .= 

= (M) < 0|T(xAf(^l)riXAf(^l) • • -XMiXn^nXMiXn) 

i 

x e 



i / d 4 x{ - gx M {x)YXM{x)C u {x) - kxm{x)Xm{x)(t{x)}, {m) 

J ) l U > conn. • l 1 - 1 ^) 



The tree cr-field approximation can be justified keeping massive magnetic monopoles very 
heavy, i.e. M 3> Mc- This corresponds to the London limit M a = 2M > Mq in the 
dual Higgs model with dual Dirac strings [11]. 

In Ref.[2] Eq. fll.13D has been applied to the computation of the magnetic monopole 
condensate < x(x)x(x)] £ > in dependence of the dual Dirac string shape represented 
by the electric string tensor £^ v {x) Eq.(l.g). The magnetic monopole condensate < 
x(x)x(x); £ > has been calculated in the tree cr-field approximation neglecting the fluc- 
tuations of the dual-vector field around the Abrikosov flux line which satisfies the 
equation 

(□ + M*)C[£(x)] = -d^(x), (1.14) 

and takes the form 

C v [£{x)} = - J dVAfi-i'.Mcja'/r^), (1.15) 
where A(x — x', Mc) is the Green function 

r d A k p~ ik ■ {x - x') 

= ! W W-v-n - (L16) 

In this paper we will calculate the magnetic monopole condensate < X{x)xi x )] £ > in 
the tree cr-field approximation but taking into account quantum fluctuations of the dual- 
vector field Cfj_ around the Abriksov flux line C u [£(x)]. An important role of such fluc- 
tuations for the formation of the interquark potential has been pointed out in Ref. [12] 
within a dual Higgs model with dual Dirac strings. 
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This paper is organized as follows: In Sect. 2 we calculate the magnetic monopole 
condensate in the tree a-field approximation and explicitly integrate out quantum fluctu- 
ations of the dual-vector field around the Abrikosov flux line. In Sect. 3 we calculate 
the contibution of the string shape fluctuations to the magnetic monopole condensate. In 
the Conclusion we discuss the obtained results. 



2 Quantum dual— vector field fluctuations 

In the tree a-field approximation we determine the magnetic monopole condensate 
<x(x)x(x);£> following [1] as 

< x(x)x(x);S >-< x(0) X (0) >= (M) < 0\t(xm(x) X m(x) 

% / d 4 z{ - gxM(z)Y 'xm{z)C v {z) - k X m(z) Xm (z)<j(z)}, m 
*e J 1 |0>y o „L (2.1) 



where the r.h.s. of Eq.( |2.1| ) should vanish at C M = o = 0. 

The time ordering operator and vacuum wave-function act on the massive magnetic 
monopole fields xm and the fields of collective excitations cr and C^. 

The calculation of vacuum expectation values of time-ordered products of the dual- 
vector fields is convenient to perform by means of the path integral method 

1 c i d 4 z C cf r\CJz)] ,,,, / 
< X (x)x(x); £>-< x(0) X (0) >=-J VC,e J h W < Q\t(x m (x) Xm (x) 



i j d A z[ - gXAi{z)YXAi{z)C u {z) - k Xm (z)xm{z)(t{z)}, (m) 



xe J I ^ ' - ' ™ '~ )|0>W, (2.2) 

where Z is a normalization factor determined as 

VC^e J . (2.3) 

The effective Lagrangian £ e fj[C M (2:)] is defined by the part of the effective Lagrangian 
Eq. flLp related to the C M -field: 

C cS [C,(z)] = ±F, u (z)F^(z) + l -MlC„{z)C»{z). (2.4) 

In order to integrate out quantum fluctuations of the dual-vector field C M around the shape 
of the Abrikosov flux line we split the C^-field into a classical field C fJi [£(z)] induced by 
the Dirac string and quantum fluctuations c M (z) around that classical background. [12]: 

C ll (z) = Cp[£{z)]+c ll (z), (2.5) 



where C lx [£(z)\ satisfies Eq.( |1.14j) , and c^z) are the fluctuations of the dual-vector field 
having a vanishing vacuum expectation value < c^(z) >= 0. Substituting the decompo- 
sition Eq.( [2.5| ) in the Lagrangian Eq.( [2.4| ) we arrive at the Lagrangian of the quantum 
fields c M (x) fluctuating around the Abrikosov flux line. 

C cS [C,(z)] = C stTing (z) + ~ c,{z) [(□ + M 2 C ) tf» - d»d v ] c u (z), (2.6) 
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where we have used Eq.( |1.14|) . The Lagrangian of the dual Dirac string jC s trmg(z) is defined 
[5,11-13] 

J d*z£ stling (z) = \M 2 C J J d'zd%E m {z)^ v {z - y, M c )£^(y), (2.7) 

where A«(* - y, M c ) = « + 2d a d l/ /M c )A(z - y; M c ). 

Since the Lagrangian Eq.( |2.5| ) is Gaussian with respect to the c^-field, we are able to 
integrate out the c^-field exactly. 

Integrating over the c^-field we reduce < x{ x )x{ x )i^ > to the form: 

< x(x) X (x);£ > - < x(0) X (0) >= (M) < 0\t(xm(x) X m(x) 
1 

x e 2 



i\ g 2 II d 4 zd 4 y [xm{z)^Xm{z)\ D^z - z', M c ) \Xu{z' )i u Xm(z> )] 



% I d 4 z{-gx M (z)YXM(z)C u [£(z)] - kxm(z)xm(z)g(z)}, 

X e J J l U >conn.) 



where D^ u (z — z',M c ) is the Green function of the free c^-field: D^ u (z — z',M c ) = 
(9tj,u + d^dv/Mc) A(z — z',M c ). Since herein we consider dual Dirac strings as classical 
objects the contribution of the Lagrangian of the dual Dirac strings Astringt^) cancels out. 

The integration over cr-field degrees of freedom in the tree approximation we perform 
following [2]. This yields 



< x(x) X (x);£ > - < x(0)x(0) >= < x(0)x(0) > ^ < 0\t( X m(x) X m(x) 

-i\g 2 ff d 4 zd A z'[ X M(z)^XM(z)}D^(z-z';M c ) [ X m(z )Yxm(z' )] 



x e 



2" 



-ig d z [xm(z)Yxm(z)] C u [£(z)\ 
X e jyj > conn , l^-yj 

For the calculation of the vacuum expectation value in the r.h.s of Eq.( |2.9| ) we assume 
that the massive magnetic monopole fields Xm{x) are almost on-mass shell. It is valid 
due to a very large mass of the monopole fields. In this case the transfer momenta are 
small compared with the mass of the dual- vector field Mq- By virtue these assumptions 
we can reduce the four-monopole interaction in Eq. fl2.9p to a point-like interaction. 

< x( x )x( x );£ > - < x(o)x(o) >= -t^ < x(o)x(o) > (M) < o\t(xm(x) Xm (x) 



AM 3 

~i I d ' z \ 7^721xm(z)Yxm(z)] 2 + g [xm(z)Y Xm (z)} C u [£(z)} 



x e 



[2M C j | >W. (2.10) 



Thus, since M 3> Mc the vacuum averaging over the massive magnetic monopole fields 
can be represented by the momentum integrals [2] related to the magnetic monopole 
diagrams depicted in Fig. 1 and Fig. 2: 



^<0\t(xm(x)xm(x) 



x e 



j [2M C J^m> (M) 
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Figure 1: Magnetic monopole diagrams describing the magnetic monopole condensate 
around the dual Dirac string without fluctuations of the dual-vector field. 



d 4 k 

tr< 



16tt 2 J 7r 2 r" \M- k + gC[£(x)} M-'k. 

1_ rd*k f 1 1 1 

~16^i ^H tT { M-k + gC[£(x)] M-k + gC[£(x)] lf11 J 

~ / ^ 2 w i \ (n_1) /-^i f i i 

x LI TO ) ( —2 ) / T2- tr i 77- r, , 77^777777^777 r. , .Ar^-M ^ 

n=l 



2M 2 / V16tt 2 ; 7 tt 2 z [M-^ + ^W] ' w M-ib 1 + 0C[£(a;)] 

n-l , 



d 4 k n , f 1 1 



ir 2 i 



1 " ' M - fc n _! + ^(x)] 7 ^" 1 M - ^ + ^[f (z)] 7 ^. 



x 77^7 / ^tr k - , 1 _ - + (2.11) 



16tt 2 J ttH {' m M-k n + gC[£(x)} 

The first term in the r.h.s. of Eq. (|2.11|) has been calculated in Ref.[2] at the neglect of 
quantum fluctuations of the dual-vector field C^, whereas the second term is fully due to 
these fluctuations. We calculate the second term keeping leading divergent contributions 
as it is accepted in the MNJL-model [1,2]. The ellipses denote the contribution of the 
diagrams depicted in Fig. 2b. This contribution is less important with respect to the 
contribution of the diagrams in Fig. 2a and below we adduce the result of the calculationof 
these diagrams without comments. The vacuum expectation value Eq.( |2.11 ) amounts to 



^<Q\t(xm(x) X m(x) 



x e 



^ z {^2\Xm{z)^Xm{z)] 2 + g [xm(z)Yxm(z)} C v [£{z)]\ 

[2M C J^|0> (M) 



conn. 



M - 2 g 2 C»[£{x)]C»[£ (x)) + ^ j - £s[Ji(M) + M 2 J 2 (M)] 



in 



X 



2 co ( i 2 ~1 ( n— 1) 

9 2 L W^ [MM) + M " h (M)] C,[£{x))C»[£{x)] 



2M Cn=1 UHO 



^2 

2 



TUT 1V± C T or < 

2Gi g 2 C„[£(x)]C»[£(x)\. (2.12) 



24vr 2 _ 2 

Here we have used the definitions of Mq given by Eq. fll.2p . 
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(a) 





(b) 

Figure 2: Magnetic monopole diagramms describing the contributions to the magnetic 
monopole condensate caused by the quantum fluctuations of the C^-field. 

Thus, integrating out explicitly quantum fluctuations of the dual-vector field 
around the Abrikosov flux line and taking into account the contribution of the scalar 
field a in the tree approximation we obtain the magnetic monopole condensate in depen- 
dence on the shape of a dual Dirac string in the form: 



< x(x)x(x);£ >=< x(0)x(0) > 
x L « 2 d <X(0)X(0)> 1V±C ^2G 1 



2 , 9 2 



M C + ~yr 2 1 ^ 

1024^ AP + — -f-^W^^^ I (2 ' 13) 



6d 



The second term in the braces of Eq. ( [2.13| ) is the result of the calculation of the diagrams 
in Fig. 2b. It may be seen that quantum fluctuations of a dual-vector field around the 
Abrikosov flux line give a substantial contribution to the magnetic monopole condensate. 
In order to retain the agreement with the results obtained within CQED [14] which testify 
the suppression of the magnetic monopole condensate in the region close to a dual Dirac 
string we have to impose the constraint Mq > g 2 /6Gi. 



S 



3 Dual Dirac string shape fluctuations 



The string shape fluctuations we define following [13,15] by X M — > X M + ^(X), where 
rj^(X) describes fluctuations around the fixed surface S swept by the shape T and obeys 
the constraint tj^(X)\qs = [13,15] at the boundary OS of the surface S. 

The magnetic monopole condensate defined by Eq. fl2.13|) and averaged over string 
shape fluctuations reads 



< x(x)x(x);£> - < x(o)x(o) > 

,2 



« 2 G! <x(0)x(0)> 



1024vr 4 



M 3 



M c + 



9 



< x(o)x(o) > 



2Gi « 1 1 



M. 



, £ 2 967T 2 



W (3.1) 



c 6Gi 

where Z s h ape is a normalization factor determined as 



shape 



i <^cS N [77] 



and <5<Sn[?7] has been calculated in Ref. [13]: 

5S^r]] = JJd 4 xd*y Va (x) O af3 (x - y) n P {y) + J d'x^x) O a (x) 
The operators o/3 (x — y) and O a (x) are given by 
0^{x-y) = ld^(x-y)SM 



2 



2 
1 
4 



<9 2 



dx a dxp 



dx a dxp 
A\(x-y;M c )£^(y) 



dx 



° 5^(x-y)£ flv (x)E^(x) ° 



dy v 



+^-g^( x )^-^ a (x)5^(x-y) 



dxp 



~ M 2 C A S, u {x) -£rA\(x - y; M c ) £ x »(y) 



dxp 



+6^(x-y)^-^(x)S tlu (x)^- 



+ ]M 2 c -^-S^(x)A xp (x-y;M c )S a %y) ° 



dy h 



(3.2) 



(3.3) 
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O a (x) 



() ^(x) + ls u Ux)J-^ a (x)-^(x)-^^(x), (3.4) 



2 dx, 



where T^ Vfi (x) is determined by 

E^(x) = \m 2 c J d A z^ x (x-z-M c )S Xv (z). 

Using Eqs.( |1.8|) and ( |1.15| ) we determine g 2 C M {?7(x)}C A1 {^(x)} as follows: 
g 2 C ll {r ] (x)}C»{ V (x)} = 



(3.5) 



g 2 Q 2 II d*a Xfl (X)d*a p ,(Y)-^-A(x-X - V (X))^-A(x-Y - V (Y)). (3.6) 



d 



d 



dx x 



OXr 



The changes of the surface elements <r Xfl (X) and da pp {Y) caused by the shifts X 
X + f](X) and Y — > Y + f](Y) have not been taken into account in the r.h.s. of Eqs.flO 
since they vanish for the straight string. Indeed, the integration over the 77-field we 
perform following [13,14] for fluctuations around the shape of the static straight string 
with the length L tracing out the rectangular surface S with the time-side T. In this case 
the electric field strength £ pv (x) does not depend on time and reads 



£(x) = e z QS(x)5(y) 



Olz + -L 



(3.7) 



where at X q = (0, 0, | L) and Xq = (0, 0, — | L) quark and antiquark are placed, respec- 
tively. Then the unit vector e z is directed along the z-axis and 6(z) is the Heaviside- 
step-function. The field strength Eq. (|3.7| ) induces the dual-vector potential 



< C(x) >= -iQ 



d 3 k k 



x e. 



4tt 3 k z M 2 +k 2 



sin (^)e^. 



(3i 



Allowing only fluctuations in the plane perpendicular to the string world-sheet, i.e. setting 
f] t (t, z) = r) z (t, z) = [13,14], we arrive at the fluctuation action S S^[r] x ,rjy] 



6Sn[tIxiVv] = 

T/2 T/2 L/2 L/2 

= J dt J dt' J dz J dz' 

-T/2 -T/2 -L/2 -L/2 



d7] ^ z) 0l {t,z\t\z') d ^^ Z,) 



dt 



dt' 



<h,M - :) 2 (t, z\t', z' ) g^g^j + Vx (t, z) O s (t, z\t>, z' ) Vx (t', zf) + (x~y) 



dz dz' 
where the operators Oi (i — 1, 2, 3) are defined by 



, (3.9) 



1 (t,z\t',z') = Q 2 



d 2 k ± 
64vr 4 



dk dk z e 



—ik (t — t') + ik z (z — z')_ 



M 2 c + k 2 z + \k\ 



M 2 -k 2 + k 2 + k 2 
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2 (M t ', 2 ') = W^ / / rffaft.e-ifatf-O +*.(»-»') c ° 2 

•/ 647T 4 7 j j\zf± - fcn 4- F 

3 (t, z\t\ z> ) = 5{t-t>) 5(z -z')Q 2 J ^± J ^ sin cos(k z z 



x j^j-g / 7 J dk dk z e-Mt - O + " *0 



M2+fc2 + fc 2 "7 647T 4 



O — oe-oo 

C* ^0 



M 2 - k 2 4- A- 2 
' ''" A (3.10) 



M2-fc 2 + fc| + k 2 



The linear terms in the 77-field expansion do not appear, since only the components S tz (x) 
and T, z t{x) survive in Eq.( j3.4|) for the static string strained along the z-axis. 

The fluctuating fields rji(t,z), where i = x,y, should obey the boundary conditions 
i]i(t,z)\ds = 0, which for the rectangular surface read [13,14] 

Vi(t,z)\ 9s = Vi (±T/2,z) = Vi (t,±L/2) = Vi (±T/2,±L/2) = 0. (3.11) 

The integration over the ^-fields should be performed with the weight 

// V Vx V Vy e l5S ^^v\ (3.12) 

•'-'shape J J 

where the measure of the integration reads 

Vr] x Vri y = Yl II dr] x (t,z)dr]y(t,z). (3.13) 

-T/2<t<T/2 -L/2<z<L/2 

Before the integration over the ^-fields we can make some simplifications of the A- 
functions. For this aim we suggest to integrate out k± keeping only the main divergent con- 
tributions as it is accepted in our effective approach [1,2]. In the region — L/2 < z < L/2 
this reduces the operators (i = 1,2, 3) to the expressions 

1 (t, z\t', z' ) = 2 (t, z\t', z' ) = %^ 5(t -t')5(z-z'), 

62lX 

3 (t, z\t',z' ) = ^ ( - ^ + ^ S(t -t')8(z-z'), (3.14) 

where Aj_ is the cut-off in the plane perpendicular to the world-sheet of the string. The 
fluctuation action becomes 

_ „ T/2 L/2 

3Q 2 A 2 r r 

5Sx[ri x ,riy] = j dt J dz[rj x (t, z) (-A)r) x (t, z) + (x *-+ y)], (3.15) 

-T/2 -L/2 

where A is the Laplace operator in 2-dimensional space-time 

d 2 d 2 

A = -at + 8? (316) 
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The common factor Q 2 A 2 L /8ir can be removed by the renormalization of the 77— fields, and 
the action of the fluctuations becomes 

T/2 L/2 



5S N [r] x ,T]y} 



dt J dz 

T/2 -L/2 



Vx{t, z)\ - -A- J 7] x (t, z) + (x^ y) 



M, 



c , 



(3.17) 



The factor 1/Mq is introduced by dimensional considerations. We have used the mass of 
the dual-vector field, since the Abrikosov flux line is localized in the region of order of 
0(1/ Mc) in the xy-plane. Of course, the final result does not depend on the parameter 
making the operator A dimensionless. 

For a static dual Dirac string and after the renormalization of the 77— fields the scalar 
product g 2 C fM {r](x)}C Jl {rj(x)} amounts to 



g 2 CMx)}C»{ V (x)}=g 2 Q 2 



J(k + q)-x V 3 Q M C A 




— 7 — n — ■ — ; sin 

47r 3 47r 3 k z q z \ 2 

7T 1 2 



sin 



1 

, 2 J M 2 +k 2 



[(K + q x )Vx(t, z) + (k y + q y )r] y (t, z)\ 



M^ + q 2 

Thus, in the static dual Dirac string approximation Eq.( |3.1| ) reads 

k 2 G, <x(0) X (0)>" 



(3.18) 



< x(x)x(x);S > - < x(0)x(0) > 

■2 



1 - 



1024vr 4 



M 3 



=< x(o)x(o) > 



c 2d k 2 1 
2 # 2 96tt 2 M 2 



M. 



c 



6G1 



rf 3 fc rf 3 g fc^. + k y q y ^ f k z L\ ^ f q z L\ 1 
47T 3 47T 3 k z q z 



sm 



sm 



2 ) M 2 +P M 2 +q* 



J(k + q) ■ x 



T/2 L/2 



X 



/ Vr] x Vr] y exp —% f dt' f dz' 

hanp J J J 



7 

" shape 

5tt 1 2 



Ml, 



Vx(t',z') 



3 QM C A 
Integrating over the //-fields we get 

< x(x)x(x);S > - < x(0)x(0) > 

,2 



-T/2 -L/2 

(** + - f ) S(z - z' ) Vx (t', z') + (x^y) 



(3.19) 



« 2 G! <x(0)x(0)> 



1024tt 4 



M 3 



2 , 



=< x(o)x(o) > 



2d « 1 



x^ 2 Q 2 



x exp 



„ , 2 r 96vr 2 M 2 

d 3 k d 3 q k ± ■ q ± ( ' k z L\ 
sm 

V 2 y 

00 L/2 



47T 3 47T 3 

.87r (fcj. + g ± ) 
' 3 Q 2 A 2 ± 



sm 



J(k + q) ■ x 



2 / M 2 +k 2 M 2 + <f 



dt' J dz'5(t - t')5(z - z' )A- 1 S(t - t')8(z -z')y (3.20) 

-L/2 ' 
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where k± ■ q± = k x q x + k y q y . 

In the integrand the Green function A~ 1 5(t — t') S(z — z') should be calculated at 
certain boundary conditions. For the open dual Dirac string the calculations should be 
performed using Dirichlet boundary conditions [13,14]. Since in this case 5{z — z' ) is given 
by 



27m 



8{z-z') = - J2 sin I —[ 



-z sin 



2nn , 



(3.21) 



the Green function A 1 S(t — t') 5(z — z' ) is defined 



O CO „ 

^5{t-t')5{z-z') = - £ J 



du e 



-iuj{t — t' ) 



2im 



2tt 2 4vr 2 n 2 
uj l - 



sin 



-z sin 



2im 



(3.22) 



L 2 



Using Eq.( |3.22| ) we reduce Eq. ( P-20| ) to the expression 

k 2 G 1 <x(0)x(0)> 



< x(x)x(x);S> - < x(0)x(0) > 
c 2Gi « 



1024tt 4 



M 3 



=< x(o)x(o) > 



M, 



# 2 96vr 2 M 2 



o 



6Gi 



x exp < — 2 



d 3 /c ePg k±- q± . ( k z L 
— n — o — ; sm 

47r 3 47r 3 fc^g* \ 2 

8tt + gl) 2 2 



Q 2 A 2 



E 



sm 



da; 



'gjA l l 

, 2 / M 2 + fc 2 M 2 + r 
1 r, ( 2irn \ 1 



J, (k + q) • x 



2tx 2 4vr 2 n 2 

oo UJ — 



sm 



(3.23) 



L 2 



By applying the Wick rotation uj — > zcj we obtain the magnetic monopole condensate in 
the form 



< x(z)x(zO;£ > - < x(o)x(o) > 

,2 



1 



/c 2 Gx < X (0)x(0)> 



1024vr 4 



M 3 



=< x(o)x(o) > 



„ r2 g 2 96vr 2 M 2 



crfc erg k±- q± ( k z L 
sm — — 



x exp 



27r 2 27r 2 k z q : 
\ l(k± + q±) 2 



sm 



q z L" 



J(k + q) ■ x 



2 J M 2 +k 2 M 2 +<f 



<p(z) 



l 2 Ai j 
where we have used the Dirac quantization condition g Q = 2n and denoted 



4 £2 

3 71 L 



E 



duj 



27m 



2ix 2 47r 2 n 2 

^ UJ + 



sm 



(3.24) 



(3.25) 



L 2 
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The function <p(z) is denned by a divergent series. Therefore, it should be regularized. The 
regularization of this function we perform in the Appendix. As it is shown the regularized 
(p(z) -function equals to zero for any z ranging the values from the interval —L/2 < zL/2. 
Thus, below we set ip(z) =0 and get 



< x(x)x(x);S > - < x(0)x(0) > 

,2 



1 



k 2 G, <x(0)x(0)> 



1024tt 4 



M 3 



=< x(o)x(o) > 



M 2 



9 



2G X «' 



„ t2 g 2 96tt 2 M 2 



X 



d 3 k d 3 q k± ■ q± _. / /c^' 
27r 2 2vr 2 fc z & 



sin 



sm 



1 



1 



J(k + q) ■ x 



2 / Mq + k 2 Mq + q- 



(3.26) 



For a sufficiently long string the main contributions to the integrals over k z and q z come 
from the momenta \k z \ ~ 2/L and |^| ~ 2/L. These values are small compared with M c 
and can be neglected in the denominators. This reduces the r.h.s. of Eq. (|3.24|) to the 
form 



< x(x)x(x);£> - < x(0)x(0) > 

,2 



1 - 



k 2 G x <x(0) X (0)> 



1024vr 4 



9 



-< x(o)x(o) > 



— — 967r2 M2 
c 6Gi 



M 3 



sm 



cos(fc 2 2;) 



c? 2 /c_l c? 2 g_L 

"2^ 2 ~^ 2 ~ (M 2 + £ 2 ) (M 2 + q 2 ) 



J (k± + q±) ■ x± 



(3.27) 



where x±_ = (x,y). Taking into account that z is in the interval —L/2 < z < L/2 we 
simplify Eq. fl3.26j ) as follows 



< x(x)x(x)]£> - < x(0)x(0) > 

,2 



k 2 G! <X(0)x(0)> 



=< x(o)x(o) > 



M 2 + ^— 2 , 
c 2d k 2 1 



Mr 



1024tt 4 M 3 

27, ft. J (k± + q±) ■ x± 



d 2 k± d 2 q± (k± ■ q±) e 



2 _ 96vr 2 M 2 77 2tt 2tt [m 2 + /c 2 ] [M 2 + g 2 ] 



(3.28) 



We can represent the r.h.s. of Eq. (|3.27 ) in the more convenient form 

« 2 G! <x(0)x(0)>" 



< x(x)x(x);£> - < x(0)x(0) > 

,2 



1 - 



1024tt 4 



M 3 



9 



< x(o)x(o) > 



MA H 2 1 

c 2Gi if 1 

„ - £ 2 96tt 2 M 2 



V 



X±. 



d 2 k ± e 1 ^'^ 
2tt M c + k 2 



(3.29) 
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where y,=? is the gradiant with respect to x±. 

Integrating over directions of the vector k± and taking the gradient we get 

k*G x <x(0)x(0)>" 



< x(x)x(x);S > - < x(o)x(0) > 

c 2d « 



1024vr 4 



M 3 



- < x(o)x(o) > 



„ , 2 g 2 96tt 2 M 2 
c 6d 



00 dkk 2 Jx(kr) 
Mq + k 2 



(3.30) 



where J\(uk) is a Bessel function and r = The integral over k can be calculated 

explicitly and reads 

r°° dkk 2 JAkr) 2M 2 C dkkJ (kr) ijr T , , 

where Ki(Mcr) is a McDonald function. 

Thus, the magnetic condensate averaged over quantum dual-vector field and string 
shape fluctuations reads 



-< x(o)x(o) > 



K 2 Gx <x(0) X (0)> 



M, 



c 



2d 



K 



Ml 



1024tt 4 



M 3 



M 2 



^ 2 96vr 2 M 2 



i^i(M c r) 



(3.32) 



6Gi 



It may be seen that due to the constraint > g 2 / QG\ the magnetic monopole condensate 
at distances close to the string r — ► becomes suppressed. For r — > the McDonald 
function K\(Mcr) behaves like Ki{Mcr) — > 1/Mcr. However, we have to emphasize that 
in such a model like the MNJL model [1,2] and a dual Higgs model [11] the region of 
distances close to the string is restricted by the constraint r > 1/A±, where Aj_ is the 
cut-off in plane perperdicular to the world-sheet of a dual Dirac string [2,5,11-13]. Due 
to Nambu [5] 1/Aj_ should be understood as a thickness of a string. Following [5,11] this 
cut-off Aj_ should be identified with the mass of the a-meson, , i.e. Aj_ = M a = 2 M. 
As has been shown in Ref. [11] this choice makes next-to-leading order corrections in 
large M a expansion to the string tension logarithmically small compared with the leading 
order contribution. Thus, the McDonald function K\{Mot) is restricted from above as 
K\(Mct) < 2M/Mq- Since the value of the condensate can be either negative or zero, 
we can impose the constraint 



M 2 + 



if 



2G t 



Ml- . 
c 6d 



g z 24tt 2 



> 0, 



(3.33) 



where we have neglected the contribution of the term of order 0(1/M 3 ). Using the relation 
k 2 = 2g 2 /3 we bring up Eq. ( |3.33| ) to the form 



M 2 > 



12tt 2 



6Gi 



(3.34) 



36tt 2 
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This relation agrees with the inequality M c > g 2 /6Gi for any g 2 /367r 2 < 1. 

At distances far from the string r —>■ oo the contribution of the string is exponentially 
suppressed as e~ 2 M c r ^ ue ^ Q ^ e Meissner effect, and the magnetic monopole condensate 
tends to the magnitude of the order parameter, i.e. < x(0)x(0) >• A similar influence of 
an electric flux tube, being an analogy to a dual Dirac string in CQED, on the magnitude 
of the magnetic monopole condensate has been observed within CQED [15]. 



4 Conclusion 

The investigation of the magnetic monopole condensate around a dual Dirac string has 
shown that the integration over quantum fluctuations of the dual-vector field C M around 
the shape of the Abrikosov flux line leads to a substantial non-positive defined contribu- 
tion. The former might change the result obtained in Ref.[2] concerning the suppression 
of the magnetic monopole condensate at distances close to the dual Dirac string. In or- 
der to retain this suppression obtained in CQED [14] we have to impose the contraint 
M c > g 2 /QG\. Since the coupling constant G\ is arbitrary, the mass of the dual-vector 
field Mq is left arbitrary to the same extent. Due to the constraint M c > g 2 /6Gi the con- 
tribution of quantum dual-vector field fluctuations to the magnetic monopole condensate 
decreases at distances far from the string, where the influence of the string is exponen- 
tially suppressed due to the Meissner effect. At infinitely large distances the magnitude 
of the magnetic monopole condensate tends to the magnitude of the order parameter, i.e., 
< x(0)x(0) >. The integration over string shape fluctuations can be performed analyti- 
cally only for the fluctuations around the shape of the static straight string of length L. 
The contribution of the string shape fluctuations smoothes the suppression of the mag- 
netic monopole condensate at distances close to the string and retains the exponential 
decrease at distances far from the string. 



Appendix. Regularization of the (p(z)— function 

The function (p(z) represented Eq.( 3.25 ) is defined by a divergent expression. Therefore, 
it is requested to regularize it. For the regularization of (fi(z) we introduce an arbitrary 
infra-red parameter \x as follows 

4 9 2 2 ^ 7 duJ 1 -2 fl™ \ ..^ 

<p(z) -»• (p(z) K = - — j X. / TT" 7-2T2 sm \—r z \- ( A1 ) 



3 7i L J 2ir 2 Air 2 n^ 2 V L 



The next step of the regularization is to apply the following integral representation: 

i^!lv 7 —p- i2nnt / L 1 1 dujdpz e%Pzt 

(A2) 



rtz) R = lim ^£ / r: e -^nt/i, / / ^ , 2 [1 - cob(2p.z)]. 

M-»0 3 7T L ~n J 2,71 J J 2tx uj + n + u 



It is easy to show that integrating over t and p z we return to Eq.(A.l). 
Then, it is convenient to decompose the integrals into two parts 

v(z) R = ^\z) R -v {2 \z) R , (A3) 
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where we have denoted 



4^2 ™ J dt_ i2imt/L 7 7 d udp z 



^ yJ 3nL±<J 2tt J J 2ix u^pl + fi 2 ' 

n ~ u -oo -oo -oo 2 

V {2) (z)k = - 9 --Y 7 -i2*nt/L 7 7 dud^ 

V [ )R 3 tt L ^ J 2vr i i 2, W 2 +^ +/1 2 oos^J 



-co — oo — oo 



vr L ~1 J 2ix J J 2n lj 2 + p 2 + u 2 



3 7r L ~i J 2ir J J 2tt ^ + pi + /i 

" — u — oo — oo — oc 

Now let us perform a summation over index n which gives 

(1) _4^2 7 dt e^l L f 7 ^ ejJ 
* [Z)R 3 tt L J Am sm(7rt/L) J J 2vr U j 2 + p 2 + l i 2 ' 1 J 

-oo v -oo -oo z r 

(% ^2 7 dt e CTt / L 7 7 do;dp z e^g + 2z ) + ej^g ~ 

V {Z ' R ~ 3 tt L J Am sin(vrt/L) J J 2n ^ + 2 + 2 ■ (^-7) 

-oo v ' 7 -oo -oc " rz r 

It is convenient to proceed to polar coordinates in the plane (co,p z ) and perform the 
integration over the azimuthal angle: 

(1) Ag 2 2 J dt e mt / L °fdppJ (pt) 



v {z)r = 3¥ I 7 IVi IM^/T) J -yTjT' (A8) 







^ 2) (4 = \^j J ^T^T l-^TtVMt + 2^)) + Jo(p(* - 2.))], (A9) 
6 ir L J Am simnrt/L) J p + u 

-oo v ' 7 ^ 

where Jo(a:) is a Bessel function. Since the Bessel functions in the integrands of Eqs.(A.8) 
and (A. 9) are even under the transformation t — > —t, the integrals become 



* {z)r = 3^lJ tJ ^T^' (A10) 

2 00 7 oo 

^ (2) {z)r = ^7 Z / 1 1 Jf? [Mvit + 2z)) + Mp{t ~ 2z))] - (An; 

-oo P P' 

The dependence of z can be removed by the shifts t + 2z — > t and t — 2z ^ t: 



Z i z 

P + /* 

-oo 



9 OO OO t/\ 

(i)M _ (2)^n _ A 9 2 /■ dt fdppj (pt) 



6 TT L J All J V + U 

-oo ^ 

As the integral over t equals to 

oo 

I dtpJ (pt)=2, (A13) 
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the functions ip^(z) R and (p^(z)^ are defined by the integral over p: 
(i)/ \ (2)/ \ Zg 2 2 7 dp g 2 2 

Substituting Eq.(A.13) in Eq.(A.3) we get 

<p(z) R = 0. (A14) 
Thus, the regularized version of the ip(z) -function vanishes for all z G [— L/2,L/2\. 
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